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In disordered systems, the localization effect of the 
conduction electron enhances the electron correlation ef- 
fects drastically, leading to conspicuous features particu- 
larly in low-dimensional systems. 1-6 - 1 The Kondo effect is 
a typical phenomenon caused by the electron correlation 
around a magnetic impurity. Using perturbative expan- 
sion, the Kondo effect in the weakly localized regime was 
investigated. 7 - ) It was shown that the Kondo logarithmic 
terms are modified into the product of new anomalous 
terms and the Kondo logarithmic ones, and that the lat- 
ters are scaled into the same Kondo temperature as that 
without randomness. 8 -* The Kondo effect in strongly dis- 
ordered systems was studied by taking account of the 
Coulomb interaction among conduction electrons. 9 -' It 
was shown that the Kondo temperature has a spatial dis- 
tribution, which leads to divergence behaviors of physical 
quantities as temperature approaches zero. Such a spa- 
tial distribution in the Kondo temperature was suggested 
from experimental results of strong broadening of the Cu 
NMR line of UCu 5 _ x Pd x . 10 ) In spite of these findings, 
the effects of strong randomness itself on the behavior of 
a magnetic impurity have not yet been fully investigated 

' from a microscopic viewpoint. 

In this paper, we study the Kondo effect in two- 
dimensional (2D) strongly disordered electron systems 

, using a finite-temperature quantum Monte Carlo (QMC) 
method. 11 ^ Let us consider the single-impurity Ander- 

, son model with on-site random potentials described by 
the Hamiltonian 

H = e *4x c ^ - t ]T c\ a c ja + e d n da 

icr <ij>a (7 

+vY,( d U 0a +H.c.) + Un dr n dl (1) 

■ a 

where random on-site potentials £j are chosen to be a 
fiat distribution in the interval [-W, W] under the con- 
dition X)i e « = , < i , j > denotes the summation of 
the nearest- neighbor sites, and n da = d^da. The system 
consists of a 41 x 41 square lattice with a magnetic im- 
purity. For W > 3.0, the conduction electron is localized 
with the localization length £(W) < 37.5. 12) Thus, the 
system for W > 3.0 is probably in the strongly local- 
ized regime in low temperature. We set the condition 
€ d + (l/2){7 = and use the parameters U = 2.0 and 
V = — 1.0 in units oft. The Kondo temperature without 
randomness can be estimated as 13 ' TS- - 0.14. 
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Fig. 1. Distribution of the Kondo temperature for W = 1.0, W = 
2.0 and 3.5. Twenty-five positions of a magnetic impurity is used 
around the center. 



Shifting the position of a magnetic impurity twenty- 
four times around the center in the same realization of 
the random potential, we calculate temperature depen- 
dence of the susceptibility of a magnetic impurity x(T) in 
< W < 3.5. Depending on the position of a magnetic 
impurity in given W, its local moment can be screened 
or unscreened by the spin of the conduction electron. 
In the former case xC^O shows a local Fermi-liquid be- 
havior, while in the latter case xCO shows power-law 
or logarithmic divergence. For W = 3.5, x(^) shows 
a local Fermi-liquid behavior at seven positions. At 
three positions x{T) shows logarithmic divergence, at 
nine positions x{T) shows weak power-law divergence, 
and at the rest six positions x(^) behaves as a free spin; 
x(T) ~ T^ 1 . The exponents of the weak power-law di- 
vergence for W = 3.5 are summarized in Table I. 

The relation between the behaviors of x{T) an d the 
local density of states for the conduction electron at the 
Fermi energy p* are summarized in Table II for W = 3.5. 
Note that p l is obtained by numerical diagonalization 
method in the case of U = V = 0. As shown in Table 
II, 0.049 < pg < 0.549 for a local Fermi-liquid behavior, 
0.105 < Pq < 0.146 for logarithmic divergence, 0.043 < 
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Tabic I. Weak divergence exponents of x(T) ~ T~ a . 



Result number #1 


#2 


#3 


#4 


#5 


#6 


#7 


#8 


#9 


a 0.952 


0.938 


0.807 


0.751 


0.597 


0.485 


0.485 


0.414 


0.343 



Po — 0.162 for weak power-law divergence, and 0.012 < 
< 0.032 for a free spin. Therefore, a smaller has a 
tendency to make the localized spin unscreened, yielding 
stronger divergence of x(T) with decrease of T. 

The results in Table II also indicate that the Kondo 
temperature Tk has a spatial distribution down to Tk = 
0, since T K = [27rx(0)] _1 . 13) We extrapolate x(0) and 
evaluate the distribution of the Kondo temperature. The 
results are summarized in Fig. 1. As W increases, the 
distribution of the Kondo temperature becomes wider 
and the weight at Tk = increases considerably. We 
have thus shown that the spatial distribution of the 
Kondo temperature can be caused only by the effects 
of a random potential. 

To obtain the local information on Kondo screening, 
we calculate the correlation function between the local 
moment and the spin of the conduction electron. When 
the local moment is unscreened, the antiferromagnetic 
correlation at the magnetic impurity is suppressed no- 
ticeably in contrast with the results for the screened 
case. 14 ) This suppression of the antiferromagnetic corre- 
lation at the magnetic impurity probably cause incom- 
plete screening and a divergence behavior of the suscep- 
tibility. 

In the disordered systems with dilute magnetic im- 
purities, where each magnetic impurity acts as a single 
magnetic impurity, the observable susceptibility may be 
obtained by averaging over the susceptibility at each po- 
sition of a magnetic impurity. We thus take an average 
over the susceptibilities at twenty-five positions around 
the center of the system. The results are shown in Fig. 
2. For W = 1.5, the average susceptibility Xav{T) shows 
logarithmic divergence in 0.05 < T < T K ~ 0.14. For 
W = 2.0 and 3.5, Xav(T') shows weak power- law diver- 
gence Xav(T) - r p—ot with a _ 0.477 (W = 2.0) and 



Table II. Relation between x(T) and p z for W = 3.5 at each 
shifted position of a magnetic impurity. The numbers in the 
top row and the column of the left-hand side represent the coor- 
dinates along the horizontal and vertical axes, respectively. The 
position (21, 21) is the center of the system. FL, log-D and FS 
denote local Fcrmi-liquid behavior, logarithmic divergence and 
behavior of a free spin, respectively. #1 — #9 correspond to the 
results for the weak power-law divergence shown in Table I. 
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Fig. 2. Average susceptibility for W = 1.5, 2.0, and 3.5. The 
broken lines are fitted by the least-squares method. For W=1.5, 
2.0, and 3.5, Xav(T) = -1.36 logT + 0.397, 1.34T-°- 477 , and 
0.867 T-°' 744 , respectively. Inset: logXav(T) versus logT. 

a ~ 0.744 (W = 3.5). Therefore, in 2D disordered elec- 
tron systems with dilute magnetic impurities, the observ- 
able susceptibility shows a non-Fermi-liquid behavior in 
low temperature. 
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